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In this work we examine the turbulence maintained in a Restricted Nonlinear (RNL) 
model of plane Couette flow. This model is a computationally efficient approximation 
of the second order statistical state dynamics (SSD) obtained by partitioning the flow 
into a streamwise averaged mean flow and perturbations about that mean, a closure 
referred to herein as the RNL^ model. The RNL model investigated here employs 
a single member of the infinite ensemble that comprises the covariance of the RNLqo 
dynamics. The RNL system has previously been shown to support self-sustaining 
turbulence with a mean flow and structural features that are consistent with DNS. 
This paper demonstrates that the RNL system’s self-sustaining turbulent state is 
supported by a small number of streamwise varying modes, which form the natural 
support for the self-sustaining process maintaining RNL turbulence. Remarkably, 
truncation of the RNL system’s support to a single streamwise varying mode can 
suffice to sustain the turbulent state. The close correspondence between RNL simu¬ 
lations and DNS that has been previously observed along with the results presented 
here suggest that the fundamental mechanisms underlying wall-turbulence can be 
analyzed using these highly simplified RNL systems. 
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I. INTRODUCTION 


The analytical intractability of the Navier Stokes (NS) equations has impeded attempts 
to develop a comprehensive understanding of the dynamics of turbulence in wall-bounded 
shear flows. This impediment to understanding resulting from the complexity of NS dynam¬ 
ics has led to extensive efforts to obtain simplifications of the NS system that still retain 
fundamental aspects of the dynamics of wall-turbulence. One useful approach for simpli¬ 
fying the NS system is to study models of reduced order. Model order reduction can be 
accomplished by Galerkin projections of the infinite dimensional NS system onto a finite 
low dimensional space. The basis functions used are generally chosen for their particular 
properties such as Fourier modes for use in Cartesian channels with diffusive dissipation, 
Proper Orthogonal Decomposition (POD) projection ^- 1 2 for economy in representing the 
structures occurring in turbulence, and balanced truncation for economical representation 
of turbulence dynamics.®® 

A second area of research aimed at understanding the dynamics underlying the NS equa¬ 
tions is the study of systems in which the complexity of the dynamics has been reduced. One 
such method is confining the turbulence to a minimal channel that reduces the complexity of 
supported perturbations and results in a simplification of the flow structures. This method 
was used to demonstrate the importance of large-scale roll and streak structures in main¬ 
taining turbulence. In particular, a minimal channel was used to show that wall-turbulence 
does not self-sustain unless the confining channel is large enough to accommodate roll and 
streak structures of sufficient spanwise and streamwise extent .® 7 

Reduction of the support of turbulence has also been attempted by seeking a skeleton 
of exact coherent structures in the phase space of the transient turbulent attractor, see 
e.g. 89 . For plane Couette flow, the first such unstable solution was the fixed point com¬ 
puted by Nagata®. Subsequent research using numerical methods has uncovered additional 
fixed points and periodic orbits for plane Couette flow . 9 11 While promising conceptually, 
the project of extending these solutions to the global turbulent dynamics has yet to be com¬ 
pleted. A related approach is to isolate and model specific flow features and interactions in 
a schematic fashion, examples of this approach include the roll/streak instability SSP model 
of Waleffc 12 ( the nine-mode truncation model of Mochlis et alP®, and the model describing 
laminar-turbulent spots and structures by Tuckerman and Barkley 14 . 
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Another method of reducing the complexity of turbulence dynamics is to simplify the 
equations themselves. One example of this approach is to study turbulence using the Lin¬ 
earized Navier Stokes (LNS) equations, which can be analyzed comprehensively using linear 
systems theory.® — The LNS equations capture a number of fundamental aspects of turbu¬ 
lence dynamics including the non-normal disturbance amplification mechanism®^, which 
has been shown to be necessary for both energy production in fully turbulent flow 25 and 
to play a key role in the bypass transition mechanism. 26 The LNS equations also provide 
specific insights into the mechanism maintaining turbulence including the role of the cou¬ 
pling between the Orr-Sommerfeld and Squire equations in generating the robust transient 
growth of streaks by the lift-up process that is an integral component of the roll/streak 
mechanism underlying the self-sustaining process (SSP) in wall-bounded shear flows.®® 
The LNS equations have also proven successful in predicting second-order statistics in these 

flows.MMi 

The current work employs a statistical state dynamics (SSD) model that incorporates a 
combination of the reduction of order and the simplification of dynamics approaches. This 
model is based on the restricted nonlinear (RNL) dynamical system, which is a second order 
closure of the dynamics of the turbulent statistical state comprising the joint evolution of the 
streamwise constant mean flow (first cumulant) and the ensemble second order perturbation 
statistics (second cumulant). The SSD is closed either by parametrizing the higher cumu- 
lants using stochastic excitation 21 37 33 or by setting the third cumulant to zero, see e.g® 9 41 . 
Restricting the Navier-Stokes (NS) equations to the first two cumulants retains the nonlinear 
interaction between the perturbation momentum fluxes and the mean flow but neglects ex¬ 
plicit calculation of the nonlinear interactions among the streamwise varying perturbations. 
Turbulence was first studied in an RNL modeling framework using the Stochastic Structural 
Stability Theory system! 42 ^, which we will refer to as the RNLoo system. In the RNLoo system 
the second cumulant is obtained by solving the time dependent Lyapunov equation for the 
perturbation covariance dynamics. This time dependent Lyapunov equation represents the 
dynamics of an infinite ensemble of realizations of the perturbation structure. As a result, 
the SSD of the RNLoo system is autonomous, which is particularly useful for obtaining an¬ 
alytical results. However, the covariance matrix dimension for a perturbation dynamics of 
0(N ) is 0(N 2 ), which limits the resolution of models that can be studied directly using 
the RNLoo model. In order to overcome this limitation in subsequent implementations of 


3 










the RNL model the covariance was estimated from a single realization of the perturbation 
dynamics! 43 * 44 * This approximation shares the dynamical restrictions of the RNLqo system 
but unlike the RNLqo model it retains small fluctuations in the perturbation covariance due 
to the perturbation covariance dynamics not being that of an infinite ensemble. It is this 
dynamical restriction that is primarily responsible for the insights into turbulence dynamics 
obtained from the RNL framework.* 43 * 44 * The existence of fluctuations in the approximate 
covariance does not greatly affect the correspondence between simulation results obtained 
using the analytically tractable RNLqo and those obtained with the numerically tractable 
RNL implementation used in this work . 43 

The goals of this work are to further investigate the dynamics of RNL turbulence and to 
examine the implications of its simplified structure for understanding wall-turbulence. The 
primary focus of our study is the mechanisms that maintain the turbulent flow and determine 
the structure of its statistical mean state. It is natural that in studying these mechanisms 
attention is focused on the dynamics of the streak and roll structures. These structures, 
which are ubiquitous in wall-turbulence, respectively comprise the large-scale streamwise 
streaks of high and low speed fluid and the associated streamwise vortices the circulations of 
which reinforce the streaks through the linear non-normal lift-up mechanism. Understanding 
the dynamical mechanisms determining the structure of the mean and perturbation fields 
and maintaining the statistical equilibrium state of wall-turbulence requires understanding 
the dynamics of these structures. The rolls and streaks in wall-turbulence are not asso¬ 
ciated with an unstable linear mode but rather, are maintained by a nonlinear instability 
process that is associated with linear non-normal growth of both the streamwise invariant 
roll and streak structures and the streamwise varying perturbation held that maintains the 
roll structure. A dynamical mechanism advanced to explain this nonlinear instability is 
referred to as a self-sustaining process (SSP). The first such process traced the origin of the 
perturbations required to sustain linear non-normal streak growth to the break-up of previ¬ 
ously generated streaks.®®! An alternative explanation of the SSP suggested that Reynolds 
stresses arising from inflectional instability of the streak maintain the roll circulations .* 7 * 49 * 5 ^ 
However, analysis of simulations subsequently revealed that most streaks are too weak to 
sustain strong infl ectional instability, which resulted in the suggestion that transient growth 
gives rise to the roll-maintaining perturbations . 48 Transiently growing perturbations have 
also been shown to maintain the roll circulation in the RNL system. However, in contrast to 
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previously proposed transient growth based mechanisms!^®], the transiently growing pertur¬ 
bations in RNL turbulence result not from random the occurrence of optimal perturbations 
associated with streak breakdown but rather from systematic parametric instability of the 
time-dependent streak.® 51 


In this work we study the RNL dynamics both by approximating the closure for the third 
cumulant of the RNLqo SSD using a stochastic forcing and by approximating the closure 
of this SSD that sets the third cumulant to zero. The main results demonstrate that RNL 
turbulence is naturally maintained solely through interactions between the streamwise mean 
flow (the k x — 0 streamwise Fourier components) and a small set of streamwise varying 
modes (i.e. the k x ^ 0 streamwise Fourier components). Moreover, we establish that a 
minimal configuration for maintaining turbulence can be obtained by limiting the streamwise 
varying wavenumber support of the RNL turbulence to a single streamwise varying mode. A 
second contribution of this work is an investigation of RNL turbulence in extended channels. 
Specifically, we show that RNL turbulence self-sustains in channels with streamwise extents 
of 96nd, where 5 is the half channel height. These results suggest that self-sustaining RNL 
turbulence continues to exist in channels with infinite streamwise extent. 

This paper is organized as follows. The next section contains a derivation of the RNL 
model from the NS equations and establishes its relation to the RNLqo system. In Section 


III we describe our numerical approach. Then in Section IV we show that the RNL system 


naturally maintains a turbulent state that is consistent with that of DNS and that RNL 


turbulence is supported by only a small set of streamwise modes. In Section [IVB| we truncate 
the streamwise varying wavenumber support of the RNL dynamics to a single mode and 
explore the interval in streamwise wavenumber over which the turbulence is sustained. We 
also explore the sensitivity of the turbulence statistics to the retained wavenumber. Finally, 
we conclude the paper and point to directions of future study. 


II. MODELING FRAMEWORK 

Consider a plane Couette flow between walls with velocities ±U W . The streamwise di¬ 
rection is x, the wall-normal direction is y, and the spanwise direction is z. Quantities are 
non-dimensionalized by the channel half-width, 5, and the wall velocity, U w . The lengths of 
the channel in the streamwise and spanwise directions are respectively L x and L z . Strearn- 
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wise averaged, spanwise averaged, and time-averaged quantities are respectively denoted 
by angled brackets, (•) = d- J^ x • dx, square brackets, [•] = j- f^ z • d z, and an over¬ 
line — t Jo * w ith T sufficiently large. The velocity field u T is decomposed into a 
streamwise averaged mean, U (y,z,t) = ( U,V,W ), and the deviation from this mean (the 
perturbation), u (x,y,z,t) = ( u,v,w ). The pressure gradient is similarly decomposed into 
a streamwise averaged mean, V P(y, z,t), and the deviation from this mean, Vp(x,p, z,t). 
The corresponding Navier Stokes (NS) equations are 

U t + U- VU +VP-^AU =-(u- Vu), (la) 

K 

Lp T U ■ Vu + u • VU + Vp — -Au = — (u • Vu — (u • Vu)) + e (lb) 

it 

V ■ U = 0, V ■ u = 0, (lc) 


where £ is a stochastic excitation used to initiate turbulence and the Reynolds number is 
defined as R — U w 8/v, with kinematic viscosity v. 

Obtaining the statistical state dynamics (SSD) corresponding to (jl| requires solving 
the infinite hierarchy of cumulant equations, see e.g. 53 ^. However, a useful and tractable 
approximation to the full SSD is obtained by closing this infinite hierarchy at second order 
by either neglecting the third cumulant or parametrizing it appropriately. We refer to 
this approximation as the RNLoo system and use the oo to indicate that the ensemble 
of realizations is infinite. Physical realizations of the RNLoo system can be obtained by 
making the ergodic assumption that the ensemble averages and streamwise averages are 
equal and consequently the first cumulant is the streamwise averaged flow. In this work we 
will approximate the second cumulant of the RNLqo dynamics by the streamwise average of 
the spatial two-point correlations obtained from a single u field realization. The resulting 
approximation to the RNLoo system is governed by the following equations: 


U, + U • VU + VP - ^AU = -(u ■ Vu), 

it 

( 2 a) 

u t + U • Vu + u • VU + Vp — -5 Au = e 

it 

( 2 b) 

v - U = 0, V - u = 0. 

( 2 c) 


We refer to equation (j2| as the RNL model. The dynamics in (J2]) correspond to parametrizing 
the perturbation-perturbation nonlinearity, u ■ Vu — (u ■ Vu) in (lb), with the stochastic 
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excitation e. These equations approximate the closure of the SSD at second order using a 
stochastic parametrization of the higher order cumulants or alternatively by neglecting the 


third cumulant, which corresponds to setting e = 0. The nonlinear equation (2a) describes 
the dynamics of a streamwise averaged mean flow driven by the divergence of the streamwise 


averaged Reynolds stresses, which we denote by, e.g. (mi), (uv). Equation (2b) describes 
the influence of the streamwise constant mean flow, U(p, z, t ), on the linearized perturbation 
dynamics. 


III. NUMERICAL APPROACH 


The numerical simulations in this paper were carried out using a spectral code based on 
the Channelflow NS equations solver. 55 56 The time integration uses a third order multistep 
semi-implicit Adams-Bashforth/backward-differentiation scheme that is detailed in 5 '. The 
discretization time step is automatically adjusted such that the Courant-Friedrichs-Lewy 
(CFL) number is kept between 0.05 and 0.2. The spatial derivatives employ Chebyshev 
polynomials in the wall-normal ( y) direction and Fourier series expansions in the streamwise 
(x) and spanwise ( z ) directions. 58 No-slip boundary conditions are employed at the walls 
and periodic boundary conditions are used in the x and z directions. Aliasing errors from 
the Fourier transforms are removed using the 3/2-rule detailed in Zang and Hussaini 59 . A 
zero pressure gradient is imposed in all simulations. Table [T] provides the dimensions of the 
computational box, the number of grid points, and the number of spectral modes for the 
DNS and RNL simulations. In both the DNS and RNL simulations we use the respective 
stochastic excitations e in ([Tj) and e in (|2]) only to initiate turbulence. In order to perform 
the RNL computations the DNS code was restricted to the dynamics of (|2]). 

We also perform a number of flow simulations where the flow dynamics of the RNL model 
are restricted to a single streamwise varying perturbation and the streamwise averaged mean 
flow. We implement these cases by introducing a damping term, £u, into the equation for 


the streamwise varying perturbations (2b) in the following manner 


Uf + U • Vu + u • VU + Vp — -Au — £u = e 

it 


( 3 ) 


The geometric parameters for this set of cases is given in Table [TT] and a description of their 
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TABLE I: Geometry for the numerical simulations. x/S, y/5 and z/5 define the 
computational domain, non-dimensionalized by the channel half-height, S. N x , N y and N z 
are the number of grid points in their respective directions. M x and M z are the number of 
Fourier modes used before dealiasing and M y is the number of Chebyshev modes used in 
each simulation. 


Case 

x/S 

y/s 

z/S 

N x 

X 

Ny 

X 

N z 

M x 

X My 

X. 

M z 

DNS-A 

[0, 4tt] 

[- 1 , 1 ] 

[0,4vr] 

128 

X 

65 

X 

128 

128 

X 

65 

X 

65 

DNS-B 

[0, 8vr] 

[- 1 , 1 ] 

[0 , 4tt] 

128 

X 

65 

X 

128 

128 

X 

65 

X 

65 

DNS-C 

[0, 12 tt] 

[- 1 , 1 ] 

[0,4vr] 

128 

X 

65 

X 

128 

128 

X 

65 

X 

65 

DNS-D 

[0, 167r] 

[- 1 , 1 ] 

[0,4vr] 

128 

X 

65 

X 

128 

128 

X 

65 

X 

65 

RNL-A 

[0, 4tt] 

[- 1 , 1 ] 

[0,4vr] 

16 

X 

65 

X 

128 

16 

X 

65 

X 

65 

RNL-B 

[0, 8 tt] 

[- 1 , 1 ] 

[0,4vr] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

RNL-C 

[0, 12 tt] 

[- 1 , 1 ] 

[0,4vr] 

48 

X 

65 

X 

128 

48 

X 

65 

X 

65 

RNL-D 

[0,16vr] 

[- 1 , 1 ] 

[0,4vr] 

64 

X 

65 

X 

128 

64 

X 

65 

X 

65 

RNL-E 

[0, 24vr] 

[- 1 , 1 ] 

[0,4vr] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

RNL-F 

[0, 32vr] 

[- 1 , 1 ] 

[0,4vr] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

RNL-G 

[0,48vr] 

[- 1 , 1 ] 

[0,4vr] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

RNL-H 

[0,72vr] 

[- 1 , 1 ] 

[0,4vr] 

138 

X 

65 

X 

128 

138 

X 

65 

X 

65 


implementation is provided in Section IV B 


IV. RESULTS 

In this section we first demonstrate the agreement of RNL simulations with DNS in 
different channel configurations. We then show that turbulence in this RNL system is 
naturally supported by a small number of streamwise varying modes. In addition, the system 
can be truncated so that it is supported by a single streamwise varying mode interacting 
with the mean flow. All simulations in this section are at R — 1000; the full parameter set 
for each simulation is given in tables [T] and |TTJ 


A. The natural support of the RNL system 

RNL system has been shown to maintain turbulence with a mean flow and structural 
features that closely resemble those of DNSJ 43 !®- 52 In the current work we study the RNL 
system described by (J 2 J across a range of channel configurations. Figure [l] demonstrates 
that RNL simulations produce mean velocity profiles consistent with those from DNS in 
channels with streamwise extents varying from AnS to 16nS. Based on the results shown 
in Figure [T| and observations that the flow structures are consistent with those reported 








TABLE II: Geometry for the numerical simulations of the truncated RNL system. Aj is 
the wavelength of the undamped streamwise varying perturbations. x/8, y/5 and z/8 
define the computational domain, non-dimensionalized by the channel half-height, 5. N x , 
N y and N z are the number of grid points in their respective directions. M x and M z are the 
number of Fourier modes used before dealiasing and M y is the number of Chebyshev 
modes used in each simulation. 


Case 

A d 

x/5 

y/5 

z/8 

N x 

X Ny 

x N z 

M x 

XMy 

x M z 

D1 

And 

[0, 

At:] 

[- 1 , 1 ] 

[0,4tt] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D2 

2n5 

[0, 

47t] 

[- 1 , 1 ] 

[0, At:] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D3 

AttS/3 

[0, 

47t] 

[- 1 , 1 ] 

[0,4tt] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D4 

7 t5 

[0, 

47t] 

[- 1 , 1 ] 

[0, At:] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D5 

4:7:5/5 

[0, 

47t] 

[- 1 , 1 ] 

[0,4tt] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D6 

2n5/3 

[0, 

47t] 

[-U] 

[0,4tt] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

D7 

47t<5/7 

[0, 

47t] 

[-1.1] 

[0,4tt] 

32 

X 

65 

X 

128 

32 

X 

65 

X 

65 

K2 

24t :5 

[0, 

487t] 

[-1.1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K3 

167 :5 

[0, 

487t] 

[-1.1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K4 

127 :5 

[0, 

487t] 

[-1.1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K5 

487t5/5 

[0, 

487t] 

[-1.1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K12 

An5 

[0, 

487t] 

[-1.1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K16 

37 :5 

[0, 

487t] 

[-1,1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K20 

12tt5/5 

[0, 

487t] 

[-1,1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K24 

2tt6 

[0, 

487t] 

[-1,1] 

[0, At:] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K30 

87t<5 / 5 

[0, 

487t] 

[-1,1] 

[0,4tt] 

96 

X 

65 

X 

128 

96 

X 

65 

X 

65 

K40 

6tt5/5 

[0, 

487r] 

[-1,1] 

[0, 4tt] 

128 

X 

65 

X 

128 

128 

X 

65 

X 

65 

K50 

24tt5/25 

[0, 

487t] 

[-1,1] 

[0, 4tt] 

160 

X 

65 

X 

128 

128 

X 

65 

X 

65 

K60 

47t<5/5 

[0, 

487t] 

[-1,1] 

[0, At:] 

160 

X 

65 

X 

128 

128 

X 

65 

X 

65 


iiP, we conclude that the RNL system captures the essential features of turbulent flow over 
a range of operating conditions. In addition, we observe that varying the channel length 
does not significantly affect the nature of turbulence sustained in the RNL system. We now 
show that RNL turbulence is naturally supported by a greatly reduced set of streamwise 


varying wavenumbers. In particular, we demonstrate that when e in equation (2b) is set 


to 0, the RNL model reduces to a minimal configuration in which only a finite number 
of streamwise varying perturbations are maintained while energy in the other streamwise 
varying perturbations decays exponentially. We refer to the set of streamwise wavenumbers 


that remain when e = 0 in equation (2b) as the natural support for the RNL system. In 


order to quantify the energy in the perturbations we define the streamwise energy density 
as the perturbation energy associated with streamwise wavelength A n , which is given by 


ExM = 


1-5 JO 


2 U a „(x,y,z,t) 2 dx dy dz. 


(4) 
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(a) (b) 

FIG. 1: Turbulent mean velocity profiles (based on a streamwise, spanwise and time 
average) from (a) simulations of the RNL model and (b) DNS for channels with 
L x = {47t, 8tt, 127r, 167 t}<5. 
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(a) (b) 

FIG. 2: Selected streamwise energy densities for cases (a) RNL-D and (b) DNS-D at R — 
1000 for a channel with L x = 16ir5. The energy densities of the streamwise varying 
perturbations that are supported in the RNL simulation are shown in the following manner 
Ai = 167r<5(black), A 2 = 87r5(red), A 3 = 167r<5/3(blue), A 4 = 47r5(cyan), A 5 = 

167 r5/5(violet), A6 = 87r<5/3(gold). The modes that decay when the RNL is in a 
self-sustaining state are shown in grey in both panels. 




Here is the perturbation, u = (u, v, w), associated with Fourier components with stream- 
wise wavelength A n . 

Figure [2] compares the natural support of the RNL system to the full streamwise energy 


spectrum supporting turbulence in DNS. Figures 2a and 2b shows the time evolution of 
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the streamwise energy densities, E\ n , for cases RNL-D and DNS-D, respectively. Both of 


these simulations were initiated with a stochastic excitation, respectively e in (lb) and e in 


(2b), containing a full range of streamwise and spanwise Fourier components. The excita¬ 


tion was terminated at t — 500. Figure 2b demonstrates that in the DNS all streamwise 
perturbations remain supported. In contrast, the number of streamwise Fourier components 
supporting RNL turbulence rapidly converges to the small number in the natural set. For 
the RNL simulation in a 167 t< 5 channel the energy densities associated with all streamwise 
wavelengths < 8/3ttS decay rapidly after the stochastic excitation is removed. These 
results demonstrate that simulations using the RNL system require substantially less com¬ 
putational resources than DNS. In what follows, we will further characterize the natural 
support of turbulence in the RNL system by examining the streamwise energy density as a 
function of the streamwise extent of the channel. 

The difference between the unforced RNL system and the NS equations is the removal 
of the perturbation-perturbation nonlinearity term, u ■ Vu — (u ■ Vu). The emergence of 
the natural set of streamwise varying modes that support self-sustaining RNL turbulence 
demonstrates the existence of a set of active perturbations that are maintained though 
interactions with the time-dependent streamwise averaged mean flow. The perturbations 
that naturally decay in the unforced RNL system cannot be sustained without excitation 
from nonlinear perturbation-perturbation interactions. The fact that the absence of these 
weakly interacting perturbations does not significantly alter either the mean flow or the 
maintenance of turbulence suggests that these perturbations do not play a significant role 
either in the maintenance or regulation of the RNL turbulent state and are in this sense 
inactive. Previous work showing the close correspondence in the mean profile and the 
time-averaged Reynolds stress components shown of Figures 1 and 6 in Thomas et al. 43 
underscores this point. 

The natural support of RNL turbulence in a channel with a streamwise extent of 167rh 
consists of all of the streamwise varying modes with wavelengths longer than 8/37 t5. This 
lower limit was previously observed to be Ati/35 for case RNL-A. Each of the RNL simulations 
described in Table [T] demonstrate a similar lower limit on the wavelengths included in their 
natural support. Based on the plot in Figure [2a] it remains unclear whether there is an upper 
bound on the wavelengths comprising the natural set. In order to investigate this question 
further we next consider the unforced RNL model in longer channels. 
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Time 


(a) 

FIG. 3: Streamwise energy densities for (a) c 
R = 1000 with streamwise channel lengths L 



Time 


(b) 

se RNL-F and (b) case RNL-G in Table [I] at 
= 327 t5 and L x = 487t< 5 respectively. 


Figures [3a] and[3b]respectively show the time evolution of the streamwise energy densities, 
E\ ri . for cases RNL-F and RNL-G with respective streamwise channel lengths L x = 32n8 
and L x = 487 t5. In both of these simulations, turbulence is initiated by applying a stochastic 


excitation e in equation (2b) from t — 0 to t=500. Figure 3a reveals that in the channel 
with L x = 327t< 5 the streamwise energy density E\ 1 is consistently lower than E\ 2 . It 
should be noted that this is the first channel length in this study where this occurs, i.e. 
in cases RNL-A, RNL-B, RNL-C and RNL-D, which are described in Table [IJ the energy 
density associated with Ai has the greatest magnitude. Figure [3b] that for case RNL-G the 
longest wavelength, Ai = 487 t 5, clearly decays to zero, confirming the existence of an upper 
wavelength limit of the natural support of RNL turbulence. The maximum wavelengths in 
the natural support of RNL turbulence in channels with streamwise extents of L x = 487t< 5 
and L x = 647t 5 are A = 247 t8 and A = 327 t8, respectively. From the data across a number of 
channel configurations we estimate that the lower and upper bounds of the natural support 
for systems at R = 1000 are approximately 2tt8 and 327 t< 5 respectively. 


Figure [4] shows the wavelength, A, at which the largest time-averaged streamwise energy 
density occurs for each of the cases described in Table [TT| For RNL simulations with L x < 
327 tS the longest wavelength has the highest mean energy density. However, for the RNL 
simulations in channels with L x > 327 t5, the most energetic wavelength is near 167 t< 5. The 
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L x //r<5 

FIG. 4: (a) The longest wavelength, X max (black squares), shortest wavelength, X min (blue 
triangles) and wavelength with the largest average energy density, A (supER™) (red circles) 
in the natural set of the RNL system versus streamwise channel length, L x , for R = 1000. 
Lines are shown solely to guide the eye. 

most prominent exception is the L x = 647rh case. In this case the energy density for Ag = 
327 t< 5/3 is 1.7% greater than the energy density associated with A 4 = I 67 r<5. We also note that 
cases RNL-E and RNL-G with respective channel lengths of 2Att5 and 72^6 do not admit 
I 67 t5. In case RNL-E the two nearest wavelengths 247 t< 5 and 127 t< 5 have nearly identical 
energy densities, and for case RNF-G the maximum energy density occurs at A 4 = 187 rh. 
We conclude that is a reasonable approximation for the ideal length for structures 

interacting with the mean flow regardless of the channel length. 

Figure [| also shows the long and short wavelength limits, respectively denoted A max and 
A m i n , of the natural set of the RNL system at R = 1000. The short wavelength limit shows 
the same relative invariance with respect to the streamwise extent of the channel, L x , as the 
wavelength of maximum energy density. For cases with L x G [ 167 T, 967t]<5 the short wavelength 
limit is constrained to a small range of values A G [57t/2, 37t]<5, whereas cases with L x < 16 ttS 
have smaller Avalues. We do not see the same invariance in the long wavelength limit of 
the natural support. For simulations with L x G [27r,247r]h the long wavelength limit, A max , 
corresponds to Ai , the wavelength associated with the channel length. For channels with L x 
G [327T, 967t]<5 the long wavelength limit, X ma x, corresponds to A 2 , the wavelength associated 
with the channel half-length. In describing the long wavelength limit, we are limited by the 
discretization in the streamwise direction. If there is an invariant for the long wavelength 
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limit, it could only be discovered by examining cases with extremely large streamwise channel 
lengths; which is a direction for future work. For simulations in channels with L x E [ 27 T , 167 t ]<5 
the largest time-averaged streamwise energy density, A(supFR n ), also occurs at wavelengths 
associated with the channel length, Ai. For systems with L x E [ 247 T , 967 t ]5 this maximum 
generally occurs at A = 167 t5 unless the streamwise discretization does not permit this 
wavelength. In those cases A(supFR n ) occurs at the permitted wavelength that is closest to 
167 t 5 , i.e. for L x = 727 t5 the maximum time-averaged streamwise energy density occurs at 
A 4 = 187t£. 


B. Truncating the RNL support 


The results in the previous section illustrate the natural support of the RNL system. 
We now show that RNL turbulence can be supported even in cases where the dynamics 
are restricted to a single streamwise varying perturbation interacting with the streamwise 
averaged mean flow. We then explore the relationship between the RNL system’s natural 
support and the support of the truncated RNL system. In Section |IV A| we observed that 


simulations of case RNL-G exhibit both the upper and lower limit of the natural set, so 
we study the truncated systems in channels with L x = 487 r 5 . The numerical details for 


the truncated RNL simulations are provided in Section |III| and a full list of the simulation 
parameters for these cases is given in Table [TT[ For all of the truncated cases the simulations 
arc initiated by applying the excitation e to ( |2b[ ) for t = 500 convective time units. We then 
apply a damping of £ in the manner shown in equation (J3]) that varies from 0 to 1 /At over 
a period of 100 convective time steps to all of the Fourier components in the perturbation 
equation (the streamwise varying k x 7 ^ 0 modes) except for the single streamwise varying 
mode that we are limiting the dynamics to. For the purpose of comparison we refer to the 
RNL system in which we do not damp any modes, i.e. £ = 0, as the baseline system, this is 
case RNL-G in Table [U 

Figure [5] demonstrates that truncated RNL systems all have an appropriately shaped 
mean velocity profile (based on a streamwise, spanwise and time average) and that these 
profiles are qualitatively similar to that obtained with the baseline RNL system. In particu¬ 
lar, for cases K3 and K4 the profiles are almost identical and have R T values of 53.4 and 57.0 
respectively, which are close to the value of R T = 56 in the baseline system. Cases K12, K16, 
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FIG. 5: Turbulent mean velocity profiles (based on a streamwise, spanwise and time 
average) obtained for cases RNL-G, K3, K4, K12, K16, K24, K40 all at R = 1000 for a 
channel with L x = 487 t 5 . Cases where the RNL is truncated to longer wavelengths exhibit 
a mean profile similar to the baseline RNL simulation (case RNL-G). 


K24 and K40 show a higher shear stress at the wall and correspondingly higher R r values 
of 64.5, 69.6 and 64.5 respectively. Selective filtering of streamwise varying perturbations in 
the RNL system was previously shown to strongly influence both the mean velocity profile 
and the spanwise spectra of the velocity field ] 60 Previous work has also shown that higher 
dissipation is associated with the inclusion of shorter wavelength perturbations . 43 This im¬ 
plies that the higher wavenumber and therefore higher dissipation cases K12, K16, K24 and 
K40 require a higher shear stress in order to attain a statistical equilibrium, which is the 
trend observed in the mean velocity profiles of Figure [5] 


Figures 


6 a 


and 


6 b 


show the time evolution of *j2E\ n for simulations of case K3 and case 
K4, respectively, in Table [TTJ The wavelengths of the untruncated modes, respectively A 3 = 
16n5 and A 4 = 12nS, are both within the natural support of the RNL model in a channel 
with L x = 487 t 5 as shown in Figures [3b] and |4| Both of these cases exhibit self-sustaining 
turbulent behavior despite the severe restriction in the streamwise harmonics. The rapid 
decay of the truncated wavenumbers subsequent to the removal of the stochastic excitation 


at t = 500 is evident in figures 6 a and 6 b 


Figures 


6 c 


and 


6 d 


respectively replot ^j2E\ n in figures 


6a 


and 


6 a 


perturbation velocity, 


B u := 


alongside the RMS 

(5) 
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(b) 




(c) (d) 

FIG. 6 : Panels (a) and (b) show yj2E\ n versus time for cases K3 and K4, respectively. 
Panels (c) and (d) show yj2E\ n for the undamped wavelength (A 3 = 167 t <5 and A 4 = 127t< 5 
respectively), the RMS perturbation velocity, E u , the RMS streak velocity, U s , and the 
RMS roll velocity, Urm for cases K3 and K4, respectively. In all panels, R = 1000. 


the RMS streak velocity, 


U s := 




(U - [l /]) 2 dy dz, 


( 6 ) 


and the RMS roll velocity, 


U Roll := A / I ~ J V 2 + W 2 dy dz. 


(7) 


We note that in both simulations K3 and K4 both the RMS perturbation velocity, E u and 
the RMS roll velocity, U ro u return to levels comparable to those prior to the truncation, 


after a transient period. From figures 6 c and 6 c it is clear that this adjustment occurs over 


16 













(a) 


1 _ 


lxio 1 r 


lxio 2 r 


lxio 3 r 


1x10" 


_ .» - ; •">-v v. *. \\ ^ 

.» 


- RMS Perturbation Velocity 

. RMS Streak Velocity 

- RMS Roll Velocity 

V 2E -l24 

I I ... I I ■ 


400 


800 


1200 


1600 


2000 


Time 

(b) 


FIG. 7: (a) ^2ER n versus time for case K24. (b) yj2E\ n for the undamped wavelength 
(A 24 ), the RMS perturbation velocity, E u , the RMS streak velocity, U s , and the RMS roll 
velocity, Urm for case K24. In all panels, R = 1000. 


a much longer time interval for case K3 as compared to case K4. 

Figure [7] shows the same results as in figure [ 6 ] for case K24. The single undamped 
streamwise varying perturbation wavelength of A 24 = 2n5 for this case is well outside of 
the natural set of the associated RNL-G (baseline) case, however the system is still able 
to sustain RNL turbulence with a reasonable mean profile as seen in Figure [5} For this 
case yj2E\ 24 is approximately 5% of the RMS perturbation velocity of the flow. When the 
damping is applied at t — 500, the value of yf2E\ 24 increases by an order of magnitude such 
that yj2E XM is nearly the same as the time-averaged RMS perturbation velocity prior to 
the application of the damping. 

Figure [ 8 a] shows the streamwise energy densities for case K2 in which only A 2 = 247t< 5 
and the mean flow are retained. In this case the flow becomes laminar after the damping is 
introduced and this behavior leads us to conclude that the RNL system in a 487 t5 channel 
cannot sustain turbulence when the dynamics are restricted to perturbation structures of 
length A 2 = 247t< 5 interacting with the mean flow. The RNL system also cannot sustain 
turbulence for case K1 in which only Ai = 487 t <5 is retained. 

As in the natural support of the RNL system, there is also a short wavelength limit 
beyond which turbulence is not robustly maintained in an RNL system that is truncated 
to a single streamwise varying mode interacting with the mean flow. This can be seen in 
figures [ 8 b] and 8 d which show a short wavelength boundary for maintaining RNL turbulence 
with a single wavenumber. As in the case of the long wavelength boundary; the undamped 
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FIG. 8 : Panels (a) and (b) show a/ 2E\ n versus time resulting for cases K2 and K60, 
respectively. Panels (c) and (d) show 2 E\ n for the undamped wavelengths (A 2 and Ago 

respectively), the RMS perturbation velocity, E u , the RMS streak velocity, U s , and the 
RMS roll velocity, Urm for case K2 and K60, respectively. In all panels, R = 1000. 


wavelength, in this case Ago, does not show any increase in energy. There is a region of 
uncertainty concerning the maintenance of RNL turbulence in cases where a relatively short 
wavelength is retained. For the conditions used in K50, K60 and K70, some cases maintained 
turbulence while others relaminarized. Similarly some of the simulations with the conditions 
of case D 6 , which corresponds to a truncated RNL simulation retaining only Ag = 27 t< 5/3, 
maintain turbulence and others relaminarize. However, in case D7 which retains only the 
mode associated with A 7 = 47 t 5 / 7 , the simulation always relaminarizes. Characterizing the 
factors that lead to the sustenance of turbulence versus relaminarization given the same 
simulation conditions is a topic of ongoing work. 

A comparison between the natural support of streamwise varying modes that support 
RNL turbulence and the single streamwise varying modes that support turbulence in the 
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FIG. 9: Time-averaged RMS perturbation velocity, E u (black squares), RMS streak 
velocity, U s (red circles), and RMS roll velocity, Urm (blue triangles) for the truncated 
RNL systems in Table [IT] for channels with (a) L x = 487r<5 and (b) L x = Air5. All 
simulations are at R = 1000. The shaded region in each graph indicates the natural 
support for the RNL model in (a) a 487 t< 5 channel and (b) a ind channel. Both panels show 
that the set of perturbation wavelengths that sustain turbulence in the truncated RNL 
system is disjoint from the natural support of the RNL dynamics. The shortest wavelength 
at which simulations of the truncated RNL sustain turbulence is much smaller than the 
short wavelength limit of the natural support of the RNL system. 


truncated RNL system is provided in figures [9a] and 9b For channels with L x = 487iA, 
the wavelengths associated with the natural support extends from 2nd to approximately 
247 t5. In contrast, the long wavelength limit associated with the truncated RNL system is 
approximately 167 tS and the short wavelength limit is well beyond the lower limit of the 
natural set and is closer to nd. A full characterization of the factors determining the bounds 
of the both the natural set of modes that support RNL turbulence and the set of single 
streamwise varying modes that can maintain turbulence through interactions with the mean 
flow will provide further insight into the nature of RNL SSP and are the subject of ongoing 
work. 


V. CONCLUSIONS AND DIRECTIONS FOR FUTURE WORK 

In this paper we have shown that the RNL system intrinsically produces a minimal rep¬ 
resentation of turbulence that is supported by a streamwise averaged mean flow and a small 
set of streamwise varying perturbations. This minimal state arises spontaneously when 
the stochastic parametrization of the nonlinear interactions among the streamwise varying 
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perturbations is eliminated in the perturbation equation. The retained modes actively par¬ 
ticipate in the transfer of energy from the time-dependent streamwise averaged mean flow to 
the perturbation field and we refer to these as the natural support of self-sustaining turbu¬ 
lence in the RNL system. We further show that the ability of the RNL system to self-sustain 
turbulence is remarkably robust. In particular, we demonstrate that RNL turbulence can 
be sustained when the dynamics are limited such that they comprise a single streamwise 
varying mode interacting with the streamwise averaged mean flow. The wavelengths of the 
streamwise varying perturbations that comprise the natural support of the RNL turbulence 
lie in a closed interval and our calculations suggest that this interval becomes independent of 
channel size for long enough channels. The set of single streamwise wavenumbers that sup¬ 
port RNL turbulence is found to extend to shorter wavelengths than those that are present 
in the natural set but not to longer ones. 

The results presented here and the previously reported close correspondence between RNL 
simulations and DNS suggest that the fundamental mechanisms underlying wall-turbulence 
can be analyzed using highly simplified RNL systems. The RNL framework provides distinct 
advantages over other minimal models that have been employed. First, the equations are 
directly derived from the NS equations and are easy to implement within an existing DNS 
code. Second, implementations of RNL minimal models do not rely on a particular Reynolds 
number or channel size and therefore Reynolds number trends as well as the dynamics of the 
RNL SSP in large channels can be explored using these models. Finally, the RNL system 
inherently captures the dynamics of key flow structures as intrinsic elements of the system 
dynamics in a computationally and analytically tractable framework. These advantages 
make it a powerful tool for probing the dynamics of wall-turbulence. The insight gained 
through such studies can then be tested using DNS and exploited to develop flow control 
strategies. 
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